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Abstract. The low energy behavior of the Adler function D{Q'^) is 
studied by employing recently derived integral representation for the 
latter. This representation embodies the nonperturbative constraints 
on D(Q^), in particular, it retains the effects due to the nonvanish- 
ing pion mass. The Adler function is calculated within the developed 
approach by making use of its perturbative approximation as the only 
additional input. The obtained result agrees with the experimental pre- 
diction for the Adler function in the entire energy range and possesses 
remarkable stability with respect to the higher loop corrections. 

1 Introduction 

The Adler function D{Q'^) [1] plays a crucial role in various issues of elementary par- 
ticle physics. Specifically, theoretical description of such strong interaction processes 
as e^e~ annihilation into hadrons and inclusive r lepton decay is inherently based on 
this function. Moreover, Adler function is essential for confronting precise experimen- 
tal measurements of some electroweak observables with their theoretical predictions. 
In turn, the latter represents a decisive test of the validity of the Standard Model and 
imposes strict restrictions on possible "new physics" beyond it. 

The perturbation theory still remains the only reliable tool for calculating the 
Adler function at high energies. Namely, in the asymptotic ultraviolet region D{Q'^) 
can be approximated by the power series in the strong running coupling as{Q'^)- How- 
ever, spurious singularities of the latter, being the artifacts of perturbative calcula- 
tions, invalidate this expansion at low energies. In turn, this significantly complicates 
theoretical description of the low-energy experimental data, and eventually forces one 
to resort to various nonperturbative approaches. 

An important source of the nonperturbative information about the hadron dy- 
namics at low energies is provided by the relevant dispersion relations. The latter, 
being based on the general principles of the local Quantum Field Theory, supply one 
with the definite analytic properties in a kinematic variable of a physical quantity at 
hand. The idea of employing this information together with perturbation theory and 
renormalization group (RG) method forms the underlying concept of the so-called 



"analytic approach" to Quantum Chromodynamics (QCD) [2J. Some of the main 
advantages of this approach are the absence of unphysical singularities and a fairly 
good higher-loop and scheme stability of outcoming results. The analytic approach 
has been successfully employed in studies of the strong running coupling [21 13], per- 
turbative series for QCD observables (see paper |4j and references therein), meson 
spectrum |5], chiral symmetry breaking [6|, and electromagnetic pion form factor [7]. 

The primary objective of this paper is to study the infrared behavior of the Adler 
function by employing recently obtained integral representation [8]. The latter has 
been derived in a general framework of the analytic approach to QCD, the effects 
due to the pion mass being retained. It is also of a particular interest to examine the 
stability of the calculated Adler function with respect to the higher loop corrections. 

The layout of the paper is as follows. In Section 2 the dispersion relation for the 
Adler function and its interrelation with the i?-ratio of electron-positron annihilation 
into hadrons are overviewed. In Section 3 a novel integral representation for -D(Q^) 
is discussed and the calculation of the Adler function within the developed approach 
is presented. In Section 4 the obtained results are summarized. 



2 The Adler function 

The Adler function D{Q^) [Ij naturally appears in the theoretical description of the 
process of electron-positron annihilation into hadrons. Specifically, the measurable 
ratio of two cross-sections is proportional to the discontinuity of the hadronic vacuum 
polarization function n(g^) across the physical cut: 

, a (e+e~ ^ hadrons: s) 1 , x 
R(s) = ^ = — lim n s - ie) - U(s + ie)] 1 

with s = > being the center-of-mass energy squared. It is worth noting here 
that R{s) vanishes identically for the energies below the two-pion threshold due to 
the kinematic restrictions, see also Ref . [9] . The mathematical implementation of the 
latter condition consists in the fact that n(g^) has the only cut > 4m^ along the 
positive semiaxis of real q'^ and satisfies the once-subtracted dispersion relation [H [H] 

ntf)^n(,S)-(,= -,S)£^j3^^<^. (2) 

where m^r — 135 MeV is the mass of the 7r° meson. 

For practical purposes it proves to be convenient to deal with the so-called Adler 
function pj, which is defined as 
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and, therefore, does not depend on the choice of subtraction point Qq in the dispersion 
relation ([2]). In Eq. ^ = — > denotes a spacehke momentum. In addition to 
the relevance to the strong interaction processes mentioned in the Introduction, the 
Adler function (JSj) plays a crucial role for the congruous analysis of hadron dynamics 
in spacelike and timelike domains. In particular, the required link between the exper- 
imentally measurable i?-ratio ([1]) and theoretically computable Adler function ([3]) is 
represented by the dispersion relation [1] 

D{Q') = Q' r T^^ds. (4) 
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At the same time, one is also able to continue an explicit theoretical expression for 
the Adler function ([3]) into timelike domain by making use of the inverse relation 

R{s) = ^ lim / Di-0 ^, (5) 



where the integration contour lies in the region of analyticity of the integrand jlOj . 

Although there are no direct measurements of the Adler function ([3]), it can be 
restored by employing the data on i?-ratio ([I]). Specifically, in the integrand of the dis- 
persion relation (jl]) one usually approximates R{s) by its experimental measurements 
at low and intermediate energies, and by its theoretical prediction at high energies. 
For the energies below the mass of the r lepton the i?-ratio (which possesses rather 
large systematic uncertainties in the infrared domain) can be substituted (up to the 



isospin breaking effects) by precise spectral function of the vector current i?^p(s) [H] 
extracted from the hadronic r decays. Thus, the i?-ratio in Eq. (j4]) can be parame- 
terized by R{s) = -R^p(s) 6^(50 — s) + -Rtheor("^) ^(^ " ^o)} where 6{x) is the Heaviside 
step-function, R^f^eori^) stands for the theoretical prediction of R{s) at three-loop 
level, 77./ = 3 is assumed, and Sq = 2.1 GeV^. The overall factor Nc Qjr is omitted 
throughout, where A^^c = 3 is the number of colors and Qf denotes the charge of the 
quark of the /th flavor. Computed in this way experimental prediction for the Adler 
function is presented in Fig. 1 by shaded band, see Ref. [12] for the details. 

As it has been mentioned above, the high-energy behavior of the Adler function ([3]) 
can be approximated by the power series in the strong running coupling as(Q^) in 
the framework of the perturbative approach. Speciflcally, at the i-loop level 

dSAQ') = 1 + [a^I\Q')y , - oo, (6) 

where ai^\Q'^) is the £-loop perturbative invariant charge. The expansion coeffi- 
cients dj are known up to the three-loop level, in particular, di = l/vr. The numeri- 
cal estimation [ISj of the uncalculated yet four-loop coefficient (i4 is adopted in what 
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follows. However, as one may infer from Fig. 1, the perturbative approximation ([6]) 
is reliable for the energies Q ^ 1.5 GeV only. Besides, expansion ([6]) is incompatible 
with the dispersion relation (jll) due to unphysical singularities of the running cou- 
pling as(Q^) in the infrared domain. The latter also causes certain difficulties for 
processing the low-energy experimental data. 



3 Novel integral representation for D(Q' 

The aforementioned integral relations (jlj) and (ISl) express the Adler function ([3]) and 
i?-ratio ([I]) in terms of each other. For practical purposes it proves to be convenient 
to express both these quantities in terms of the common spectral function. This 
objective can be achieved by employing Eqs. (jl]) and ([5]), the parton model prediction 
Ro{s) = 9{s — 4m^) [9j, and the fact that the strong correction to the Adler function 
vanishes in the asymptotic ultraviolet limit Q'^ oo. Eventually one arrives at (see 
Refs. P [12] for the details) 
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Here the spectral function Pd(o") can be determined either as the discontinuity of 
the theoretical expression for the Adler function across the physical cut Pd(c") = 
ImZ}(— cr + i0+)/7r or as the numerical derivative of the experimental data on i?-ratio 
Pd((t) = —d R{a) /din a. It is worth noting that Eq. ([7]) embodies the nonperturbative 
constraints on Adler function arising from the dispersion relation (jlj). Besides, Eq. ([8]) 
by construction properly accounts for the effects due to the analytic continuation of 
spacelike theoretical results into timelike domain. 

In order to compute the Adler function in the framework of the approach at hand, 
one first has to determine the spectral function Po{<^)- In what follows we restrict 
ourselves to the study of only perturbative contributions to the latter, namely 
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where Dp^^^^{Q'^) is given by Eq. ([6]). It is worthwhile to mention also that in the limit 
of massless pion (m^ = 0) the obtained expressions for the Adler function ([7]) and 
R-TdXio ([8]) become identical to those of the so-called Analytic perturbation theory 
(APT) [4J, the definition being assumed. 

For the illustration of the significance of the pion mass within the approach at 
hand, it is worth presenting the Adler function ([7]) computed by making use of the 
spectral function iQ for both, massless and massive cases. The obtained results 
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Figure 1: The Adler function ^ (solid curves) calculated by making use of the spectral 
function ([9]) in the massless (A) and massive (B) cases. Numerical labels correspond to the 
loop level considered. The experimental prediction for D{Q^) is shown by the shaded band, 
whereas its perturbative approximation is denoted by the dashed curve. 

are presented in Fig. 1 by solid curves. In the case of the massless pion (which is 
identical to the APT [3]), one arrives at the result, which is free of infrared unphysical 
singularities, but fails to describe the Adler function for the energies Q < 1.0 GeV, 
see Fig. 1 A. It is worth noting here that in the framework of the massless APT the 
infrared behavior of D{Q'^) can be further improved by additionally invoking into 
consideration relativistic quark mass threshold resummation [H] or vector meson 
dominance assumption [15j. 

At the same time, as one may infer from Fig. 1 B, for the case of the nonvanishing 
pion mass the representation ([7]) is capable of providing an output for the Adler 
function, which agrees with its experimental prediction in the entire energy range [8]. 
Moreover, the Adler function ([7]) is remarkably stable with respect to the higher 
loop corrections. Namely, the relative difference between the £-loop and (£ + l)-loop 
expressions for D{Q'^) is less than 4.9%, 1.5%, and 0.3% for £ = 1, £ = 2, and 
£ = 3, respectively, for < < oo, see Ref. [I2] for the details. It is worthwhile 
to mention also that the obtained results are supported by recent studies of meson 
spectrum in the framework of the Bethe-Salpeter formalism [5]. 

4 Summary 

The infrared behavior of the Adler function is examined by employing representa- 
tion ([7]), which accounts for the pion mass effects. The approach at hand possesses 
all the appealing features of the massless APT Namely, it supplies a self-consistent 
analysis of spacelike and timelike experimental data; additional parameters are not 
introduced into the theory; the outcoming results possess no unphysical singulari- 



5 



ties and display enhanced higher loop stability. In addition, the developed approach 
provides a reliable description of the Adler function in the entire energy range. 

Partial support of grants RFBR 05-01-00992 and NS-5362.2006.2 is acknowledged. 
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